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A simple symmetry analysis of in-plane and out-of-plane transport in a family of high temperature 
superconductors is presented. It is shown that generalized scaling relations exist between the low 
frequency electronic Raman response and the low frequency in-plane and out-of-plane conductivities 
in both the normal and superconducting states of the cuprates. Specifically, for both the normal and 
superconducting state, the temperature dependence of the low frequency Big Raman slope scales 
with the c— axis conductivity, while the Raman slope scales with the in-plane conductivity. 
Comparison with experiments in the normal state of Bi-2212 and Y-123 imply that the nodal 
transport is largely doping independent and metallic, while transport near the BZ axes is governed 
by a quantum critical point near doping p ~ 0.22 holes per Cu02 plaquette. Important differences 
for La-214 are discussed. It is also shown that the c— axis conductivity rise for T ^ Tc is a 
consequence of partial conservation of in-plane momentum for out-of-plane transport. 

PACS numbers: 74.25.Jb, 71.27.+a, 78.30-j 



I. INTRODUCTION 

The strong anisotropy of in-plane {ab) and out-of-plane 
(c) transport in the cuprate systems revealed by angle- 
resolved photoemission spectroscopy (ARPES), NMR, 
resistivity, Hall, Raman, and optical conductivity mea- 
surements is as unresolved and longstanding a problem 
as superconductivity itselSi*2i2iii^iS4L&. As a function of 
hole doping per Cu02 plaquette p the afc-plane resistivity 
Pab{T) (Fig. lA) shows a metallic temperature depen- 
dence {dp/dt > 0) for a wide range of doping while the 
c-axis resistivity Pc{T) (Fig. IB) varies as T"^ with an ex- 
ponent r that changes from 2 to -2 as p decreases. The re- 
sistivity ratio pc(T) / pab{T) is large and becomes increas- 
ingly temperature dependent in all (hole-doped) cuprate 
systems for p below ~ 0.22 at low temperature3242i4iSiSiL&. 

It was pointed out early on that the c— axis prop- 
erties provided an useful spectral tool to examine 
in-plane charge dynamics^. As a result, many ap- 
proaches have been put forward to address the na- 
ture of in-plane versus out-of-plane transport in terms 
of anisotropy of the in-plane quasiparticle (qp) self 
energies S](k, T), c-axis hopping ii.(k), impurity as- 
sisted hopping, interband transitions, or deconfinement 
of electrona9ii°i"ii^ii-^ii4ii5,i6,i7,i8,i9,20,2i^ Recently the 

issue of spectral weight transfers in optical conductivity 
measurements brought about by superconductivity has 
attracted a great deal of attentionZi2i2ii£. The mecha- 
nism by which 3D superconducting phase coherence sets 
in is of continued interest and debate which has been 
guided in a large part by the measurements of the c— axis 
transport properties. 

The issue is still largely unsettled basically due to the 
open question of whether electron hopping in the out-of- 
plane direction is coherentSiSiSiiSiii. If there were an at 
least partial conservation of the in-plane momentum for 
qp tunnelling along the c-axis, LDA^^ would indeed pre- 
dict an interrelation between c-axis transport and the qp 
scattering rate close to (7r,0) in the Brillouin zone (BZ). 
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FIG. 1: Experimental results for Bi-2212 for pabiT) (Panel 
A), pciT) (Panel B), the Raman-derived B2g,B\g qp relax- 
ation rate V2,i (Panel C, Panel D), respectively. The solid 
lines, circles correspond to underdoped samples (p = 0.10) 
with Tc ~ 57K, dotted lines, squares correspond to optimally 
doped samples (p = 0.15) with Tc ~ 92K, dashed lines, di- 
amonds correspond to slightly overdoped samples {p = 0.19) 
with Tc ~ 82K, and the dotted-dashed lines, triangles corre- 
spond to overdoped samples {p — 0.23) with Tc ~ 52K. All 
resistivities were measured in Ref.^, except for the overdoped 
(Tc = 52K) sample which was measured in Refi-. The Raman 
qp relaxation rates are taken from Refi— . 



What would be extremely useful would be a transport 
measurement beside conductivity which might directly 
test whether transport in the plane is intimately tied to 
out-of-plane transport. 

A behavior similar to the resistivity anisotropy is 
reflected in electronic Raman scattering measurements 
when comparing the temperature dependence of the low 
energy continuum measured in Big polarization orien- 
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tations, which project out charge fluctuations near the 
BZ axes, to configurations, which probe charge fluc- 
tuations along the BZ diagonals. Hackl et al^ and 
Blumberg and Klein^^ have pointed out the close con- 
nection between B2g Raman and in the ab-plane con- 
ductivity. Further, Opel et al?^ and Venturini et al?^ 
compared the Raman relaxation rate in each channel, 
defined as the inverse of the slope of the low energy Ra- 
man response = limQ^Q[dx'^ ^{^,T) / dVl]~^ . For 
both YBazCuaOrLi (Y-123) and Bi'sSrsCaCuzOg+a (Bi- 
2212), it was found that for i?2g symmetry, F^^ (Fig- 
IC) approximately scales with pab{T) over a wide dop- 
ing range, while for Big, Ff (Fig. ID) was found to 
cross over from metallic to insulator behavior for p less 
than ^ 0.22, occurring a higher dopings than that usu- 
ally attributed to the formation of a pseudogap^. This 
has recently been interpreted as evidence for an under- 
lying quantum critical point lying near pc — 0.22 of an 
unconventional metal- insulator transition (MIT)^^. 

At low frequencies for underdoped systems, <Jc{T) for 
Y-123 and YBa2Cu408 (Y-124) decreases rapidly with 
decreasing temperature^. From this a pseudogap has 
been inferred and well-documented. A much weaker 
spectral weight reduction is seen for La2_a;Sra;Cu04 (La- 
214)^^. Contrary to the out-of-plane conductivity, it is 
widely believed that there is no direct indication of a 
pseudogap in (Jai^- The weak dependence with temper- 
ature of the a5-plane optical sum rule compared the rapid 
decrease at low temperature of the integrated c— axis con- 
ductivity has been interpreted in Ref.i4 qps located near 
(tt, 0) from the anisotropy of t^. 

Raman scattering has been widely used to address the 
pseudogap. Recently, the presence of a pseudogap has 
been derived from c-axis Aig Raman measurements in 
Y-1242&. A much weaker signature of a pseudogap is 
seen in the Big channel in Y-123 and Bi-22132^. In 
optimal and overdoped systems, pair-breaking features 
appear only when superconducting coherence is estab- 
lished and their location at different energies for different 
symmetry channels has been well documented and in- 
terpreted as Cooper pairs having (1^2 _y2 symmetry and 
well defined low energy qp excitations^SiS^. While the 
Big pair-breaking feature appears at and scales with 
Tc for all dopings considered, closer to optimal dop- 
ing and for underdoped systems, low frequency spectral 
weight is lost at low temperatures and the Big pair- 
breaking peak becomes difficult to distinguish from the 
background^5,30,3i,32,33,34,35,36,37,38^ This loss of Spectral 

weight with temperature is very similar to the behav- 
ior seen in Raman scattering in Kondo and mixed- valent 
insulators and is indicative of gapped excitations^^. 

In the superconducting state, the temperature depen- 
dence of the a6-plane low frequency (or regular part of the 
dc) conductivity^Mi typically shows a peak around 35K 
which is material dependent and has been attributed to 
the rapid collapse of the quasiparticle (qp) inelastic scat- 
tering rate below Tc and the rise of the qp elastic scat- 
tering rate for low 7;^242iiiii^. A similar peak seen in 



in-plane thermal conductivity measurements was found 
to be sensitive to annealing conditionsi^. The c— axis 
low frequency conductivity in YBa2Cu306.95 meanwhile 
does not show a peak in this region but has an upturn 
at temperatures below 25K. The origin of the peak is 
currently not understoodi^*^. The c— axis thermal con- 
ductivity was found to show a very weak peak also sensi- 
tive to annealing conditions'^. Much less is known about 
the temperature dependence of Raman scattering in the 
static limit in the superconducting state, although some 
theoretical treatments have appeared^M^. One would 
like to test whether features shown in conductivity mea- 
surements are found in Raman scattering measurements 
and vice-versa. 

Shastry and Shraiman have noted the close similar- 
ity between the conductivity and the Raman response 
and have suggested a scaling relation exists between the 
two which follow the same temperature and frequency 
dependence^: 

n<j'{n,T)^Ax"{n,T), (i) 

with A a constant independent of frequency and tempera- 
ture. This Shastry-Shraiman (SS) relation holds if the qp 
self energy S is independent of k and has been shown to 
be exact for both the Falicov-Kimball^ and Hubbard^ 
models in the limit of large dimensions where the self en- 
ergy and vertex corrections are local. Generally though 
any k-dependence of E and/or the irreducible Raman or 
current vertices invalidates the SS scaling relation mak- 
ing it inappropriate for strongly anisotropic systems such 
as the cuprates. 

However, an approximate scaling relation may hold for 
certain cases and one purpose of this paper is to point 
out some of the connections between the conductivity and 
Raman response for strongly anisotropic systems and de- 
rive appropriate scaling relations. In particular we will, 
based on symmetry arguments, determine that a vari- 
ant of the SS relation can be formulated to show that a 
scaling relation exists separately between aab and 1/F^ 
and between Cc and as a consequence of the mo- 

mentum dependence of i±(k), in-plane self energy S(k), 
and a d^^-y^ energy gap A(k). These scaling relations 
are also found to also hold in the superconducting state. 
Comparison with the available data on Y-123 and B-2212 
in the normal state suggests that qps located near the 
BZ axes or "hot spots" become gapped above optimal 
doping^-S while the qps located along the BZ diagonals or 
"cold spots" are largely doping independent and remain 
metallic. Thus the c-axis transport is partially infiuenced 
by a correlation gap near (tt, 0) because of partial conser- 
vation of the in-plane momentum in c-axis transport and 
not completely by c-axis diffusion. There are important 
differences however with La-214. Various models for qp 
scattering as a function of doping are discussed, and it 
is found that generally no single model can adequately 
capture the complex nature of electron dynamics over a 
wide range of doping. Features of the theory in the su- 
perconducting state qualitatively describe the behavior 
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seen in the c— axis conductivity, but there are important 
questions left unanswered. In conclusion, experimental 
evidence in both the normal and superconducting states 
suggest that the in-plane momentum is at least partially 
conserved in c-axis transport over a very wide doping 
range. 

The plan of the paper is as follows: Sections II and 
III present the formalism used and the results for the 
temperature dependence of the low frequency in-plane 
and out-of-plane conductivity and Raman response in 
the normal and superconducting states, respectively, for 
the common model where t±^ vanishes along the BZ di- 
agonals, summarized in the Appendix. The results are 
summarized and open points are discussed in Section IV. 



II. NORMAL STATE 
A. Formalism 

The quantum chemistry of the tetragonal Cu-0 sys- 
tem yields an out-of-plane hopping which is modulated 
by the in-plane momentum in such a way that it is 
strongly governed by qps located along the BZ axes 
as opposed to qps along the zone diagonal, t±(k) = 
t°j_[cos{kxa) — cos(fcya)]^, as reviewed in the Appendix. 
This form for the hopping has been widely used to study 
the penetration depth^"^, c-axis conductivitjii^ii^, and 
bi-layer splittingS^ in ARPE^. However we note that 
inclusion of the Cu-0 chains or O displacements would 
lower the symmetry with the consequence that the out- 
of-plane hopping would no longer vanish along the BZ di- 
agonals which could only be noticeable at very low tem- 
peratures. We now explore the consequences of such a 
term on the regular part of the dc conductivities and the 
symmetry-dependent electronic Raman response for qp 
scattering in the a — b plane and along the c— axis in the 
following sections. 

In linear response theory, expressions for the regu- 
lar part of the conductivity and Raman response in 
the absence of vertex corrections are given as (here and 
throughout we set /cg = ?i = 1) 

Here / is the Fermi function, G^'^ are the retarded, ad- 
vanced Green's functions, respectively, = is the 
current vertex for direction a given in terms of tfie band 
dispersion and electron charge e, and is the Raman 
vertex set by choosing the incoming and outgoing light 
polarization vectors. 

The inclusion of vertex corrections is crucial for sat- 
isfying Ward identities for the conductivity and parti- 
cle number conservation for the charge density response. 



They convert scattering lifetimes into transport lifetimes, 
and also add an additional source of momentum and tem- 
perature dependence to the corresponding response func- 
tions. Vertex corrections have been recently been consid- 
ered in FLEX treatments of the Hubbard model^^ and a 
spin-fermion model^ where it was shown the Big Ra- 
man irreducible vertex is highly renormalized near the 
(tt, 0) regions of the BZ. In addition vertex corrections 
have been calculated exactly in the limit of large dimen- 
sions for the Falicov-Kimball model, where it was shown 
they are important in the Aig channel to properly lead 
to gauge invariance and particle- number conservation but 
do not contribute to other channelsi^. Generally, vertex 
corrections have not yet been generically or systemati- 
cally investigated in 2D and we thus neglect them since 
we are interested in exploring simple symmetry proper- 
ties of the various experimental probes. 

The current vertices are simply = Vx sin(fca;o), and 
ik ~ Vz[cos{kxa) — cos(fcya)]^, where Vx ^ t and Vz ''-^ 
have only a mild momentum dependence. In the limit 
where the incident and scattered photon energies are 
small compared to the bandwidth the Raman vertex is 

given as the curvature of the band: 'ja.fi — ^ '^^fc^ ™- The 
vertices are thus determined from the above band struc- 
ture as 7k — 6i [cos(A:a;a) — cos{kya)],b2sin{kxa) shi{kya) 
for Big,B2g orientations, respectively, while for c~axis 
Aig Raman 7k = a^z cos{kzc)[cos{kxa) — cos{kya)]'^ . The 
prefactors hi i, 62 ~ t' , and Qzz ~ can also be 
assumed to be only mildly frequency dependent corre- 
sponding to off-resonant scattering and therefore are only 
multiplicative constants. Since the energy range consid- 
ered is very small in comparison to all electronic band- 
widths involved the assumption 61^2 and azz to be con- 
stant is robust under all realistic circumstances. 

As can be seen by the weighting of the vertices, we 
may expect similar behavior for the B2g Raman and in- 
plane conductivity, and the -Big, c— axis Aig Raman, and 
the out-of-plane conductivity as well. The former two 
quantities assign weight around the Fermi surface (FS) 
to the diagonals while the latter three assign weight along 
the zone axes. 

In correlated electron systems the density of states 
(DOS) plays a strong role in determine transport prop- 
erties. In Mott insulators, charge transport occurs via 
excitations across a Mott gap from the lower to upper 
Hubbard bands, while in metallic systems the DOS near 
the Fermi level plays the dominant role in low frequency 
transport. The nature of how the density states evolves 
across a MIT has been an issue of intense debate for 
a large number of years as few exact results are avail- 
able. However, in the limit of large dimensions dynami- 
cal mean field theory has a great deal of insight for some 
model Hamiltoniang'^'^ . Away from half- filling the Hub- 
bard model and the Falicov-Kimball both possess metal- 
lic ground states. The DOS has a typical three-peak 
structure: the separated upper and lower Hubbard bands 
and a qp DOS at the Fermi level emerging from the 
Abrikosov-Suhl resonance in the related impurity prob- 
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lem. As the system approaches half-filhng and/or for 
larger values of U at fixed filling, the qp DOS generally 
diminishes and vanishes in the Mott insulating phase as 
spectral weight is transferred into the Hubbard bands. 
Capturing this transfer in models in realistic dimensions 
is one of the most important and difficult problem in 
condensed matter physics. 

We thus consider charge transport in correlated sys- 
tems having coherent qps as well as large energy inco- 
herent charge excitations related to the Hubbard bands. 
We model coherent qps near the FS by a phenomcnolog- 
ical momentum, frequency, and temperature dependent 
self energy derivable in principle from a renormalizable 
effective low energy theory: G^'^j^(w) = Zk(w,T)/(tj — 
Ek ± iTk(uj,T)). Here e is the renormalized band struc- 
ture, Zk(w, T) = [1 - dT,[^{cj, T)/duj]-^ is the qp residue, 
and Tk(Lu, T) is the momentum, frequency, and tempera- 
ture dependent qp scattering rate. The full Green's func- 
tion also includes an incoherent part Gmc accounting for 
larger energy excitations such as those involving the lower 
and upper Hubbard bands. In what follows we focus on 
low frequency transport in metallic phases and neglect 
Ginc and singularities of the self energy indicative of an 
incipient phase transition. 

Converting the momentum sum to an integral over an 
infinite band we obtain in the limit of low frequencies 



0,T)/dn j 



-2Nf / dx 



dfjx) 
dx 




(3) 



where Np is the density of states per spin at the Fermi 
level and (• • •) denotes performing an average over the 
FS. It can be immediately seen that the SS relation Eq. 
(|H) follows if F is independent of momentum, as it is in 
local theories^2iSiiS. In what follows we neglect specific 
features on and off the FS (such as van Hove) and ap- 
proximate the 2D FS as a circle and expand the c— axis 
dispersion for small t']^ to obtain: 

XX Conductivity, j'^ = vp sin(0), 
zz Conductivity, = v^ cos^(20). 
Big Raman, 73^^ = bi cos(20), 
i?2g Raman, 733^ = b2 sin(20), 
zz Aig Raman, 7Aig = sl^z cos^ {2(/)) (4) 

We note that the c— axis conductivity and dx'^_^ /dft 
are given by the same expressions, in accordance with the 
qp scattering rate not having a dependence. Therefore 
we confirm the SS relation for the c-axis Aig Raman and 
c-axis conductivity, respectively: 



(5) 



At low temperatures we find from Eq. ^ 



dxl^iT)/dn 




(6) 



showing the interplay of anisotropics of the scattering 
rate and the vertices governing the response functions. 

The simple expressions for a and dx" /d^ allow for 
a straightforward comparison of models for the qp scat- 
tering rate. We choose a generic model which describes 
strong scattering weighted largely along the BZ axes plus 
a temperature dependent scattering rate taken to be uni- 
form around the FS: 



rk(r) = r;,(r)cos2(20) + r,(r). 



(7) 



This form for the qp scattering rate has been widely em- 
ployed in a number of models differing in the represen- 
tations of Th and Fc constrained only to possess the full 
symmetry of the lattice (Aig)iSiiLi2iiS. Further param- 
eterizations of the anisotropy do not lead to appreciable 
differences. For the i?2g Raman as well as the in-plane 
conductivity, the vertices weight out regions of the FS 
where the scattering rate is small, along the FS diagonals 
or "cold spots". However, the Big and c— axis Aig Ra- 
man and the out-of-plane conductivity assign no weight 
to the diagonals and thus will be governed by the scat- 
tering at the "hot spots" . 

Neglecting the k-dependence of the qp residue Zk = Z, 
the resulting integrals can be easily performed to give 



on 
an 



1 



^2F,(T) ^i + F,(r)/re(r)' 

1 \ 1 



(8) 




v/i + r„(r)/r,(T) 



(9) 



= h. 



Vl + r,,(T)/F,(T)^ 
2 NpZ'^ 1 



(10) 



'2r;,(T) ^i + F,(r)/r,(r) 
|i- Vi + rh(r)/re(r) + 



F^(T)/F,(r) 



(11) 



independent of the form for Fk. 



These results for the ab-plane and c-axis conductiv- 
ity have been derived several times, most recently by 
Refsii^ii^. However here it can be seen that there is a 
direct connection between conductivities and Raman re- 
sponse functions. It is clear that the function form for the 
scattering rate determines the temperature dependence 
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of all four response functions, and that the SS relation 
Eq. does not hold in general. 

Early on, ARPES measurements yielded Tc ^ Th from 
smeared spectral functions seen near the BZ axes com- 
pared to the BZ diagonals^. However, recent ARPES 
measurements indicated bi-layer splitting may have led 
to an overestimation of Tfr^, but still the limit Fc ^ Th 
is a useful limit to explore. In this limit the response 
functions are 



2 NpZ^ 



2rvr)' 



(12) 
(13) 
(14) 
(15) 



This directly shows the similarity between the Big Ra- 
man slope and the c— axis conductivity, and B2g Ra- 
man slope and the in-plane conductivity, regardless of 
the functional form chosen for two contributions to the 
qp scattering rate. Thus in this model consistent with 
experiments, a variant of the SS relation for the cuprates 
may be expressed as 



(16) 



This demonstrates how out-of-plane transport can be di- 
rectly inferred from in-plane optical transport measure- 
ments. Further, this confirms the behavior shown in Fig. 
([T]). indicating that the in-plane momentum must be at 
least partially conserved for transport perpendicular to 
the Cu02 planes. Moreover, with Eq. this indicates 
that the c-axis Aig Raman should scale with Big Raman: 



lim x"a. 



(17) 



Eqs. H16I17|I are the central results of this section. 

When and how might the scaling relations Eqs. I|16I17|I 
breakdown? Clearly these scaling relations result from 
the momentum dependence of the respective response 
vertices, and since they are dictated solely on symmetry 
grounds, changes in how one represents the momentum 
dependence of the vertices can only lead to qualitative ef- 
fects. However, there are a number of important factors 
to consider. First, the inclusion of Ginc will change the 
scaling relations if there is appreciable spectral weight 
near the FS, but if we restrict ourselves to metallic sys- 
tems and low frequencies, then these changes are ex- 
pected to be small. They might however be large for a 
system lying near a quantum critical point and the scal- 
ing relations may be violated. Next, relating the c-axis 
conductivity to the Aig c-axis and Big Raman requires 



that the c-axis coherent hopping vanishes along the BZ 
diagonals. Deviations would come from incoherent dif- 
fusive hopping, or more complex coherent hopping paths 
such as via the Cu-0 chains in Y-123, and would result 
in a mixing in the scaling properties for in-plane conduc- 
tivity and B2g Raman transport. Lastly, vertex correc- 
tions can appreciably alter the scaling relations. Vertex 
corrections yield transport scattering rates in place of 
2E" as required by Ward identities, and yield a /—sum 
rule for the integrated conductivity proportional to mi- 
nus the kinetic energy. Ward identities can be useful 
for the conductivity to show that vertex corrections van- 
ish for a momentum-dependent self energy, but no Ward 
identities exist for Raman with crossed polarization vec- 
tors. For example, vertex corrections may renormalize 
even-parity momentum charge vertices (Raman) but not 
odd-parity current vertices (conductivity). If these scal- 
ing relations are found to hold, they would imply that 
vertex corrections at low frequencies and c-axis hopping 
along the BZ diagonals would play only a very minor role 
in determining low frequency transport. 



B. Transport models 

The scaling relations of Eq. H16|l can be seen from Fig. 
n to be qualitatively obeyed. We now consider several 
models for Tfi{T) and TdT) to explore the scaling re- 
lations Eas. H12ll5|) to address the role of anisotropic qp 
scattering. In all models, r^(T) and Tc{T) are generally 
constrained by the estimated width of the spectral func- 
tion measured in ARPES experiments^. In both a "cold 
spot"^^ and "hot spot" model^^, TdT) describes weakly 
renormalized qp scattering primarily along the FS diag- 
onals generally of the form 



-TVTo, 



(18) 



where Timp represents elastic impurity scattering and Tq 
is the energy scale of a renormalized Fermi liquid. The 
impurity scattering may be chosen to reproduce the ex- 
trapolated T — resistivity and Tq is a parameter to be 
chosen to fit a cross-over from to T in the resistivity. 
In the "hot spot" modeliiS, ThiT) = VT\7r represents 
scattering with exchange of antiferromagnetic reciprocal 
lattice momentum Q which has been widely employed 
to determine the optical conductivity, in-plane and Hall 
resistivity in relation to ARPES. However similar behav- 
ior is also obtained for scattering in systems lying near a 
charge ordering instability^ or near a FS Pomeranchuk 
instability^^. In the "cold spot" modeliSi, Th{T) is taken 
to be a constant Ths presumed to arise from strong d^-i^yi 
pairing fluctuations, and has been employed in several 
works to describe in-plane and out-of-plane optical con- 
ductivity and magneto-transporl!i2iiSiiS4i&. However the 
microscopic origin of Ths is unclear in this model. In 
the marginal Fermi liquid (MFL) model most recently 
described in RefA Tc{T) ~ T and Th{T) - constant 
due to impurity scattering in correlated systems whereby 
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TABLE I: Summary of the low temperature dependence of the 
inverse conductivities, the Raman relaxation rates and the 
scattering ratio defined in the text. 



Response 


MFL 


"Ifot spot" 


"Cold spot" 


T^{T),o-^{T) 


yl/2 


y5/4 


T 


rf(T),a-i(r) 


constant 


rpl/2 


constant 




y-l/2 


y~3/4 





strong correlation nearby a point-like scatterer induce 
real-space extensions of the impurity potential^*^. 

Following RefiSSi, the "Raman scattering rate" r^(r) 
for each channel is defined as the inverse of the Raman 

■9x''(fi-*0,T)- 



slope r^(T) = "^^^go "'" ' ^'^ order to obtain infor- 
mation on the single particle scattering rate on regions 
of the FS selected by polarization orientations /i = 1,2 
for Big 2g, respectively. In the "hot spot model" we ob- 
tain Ff ~ and F^ - T'^/'', respectively, while in the 
"cold spot" model we obtain Ff constant and F^ ^ T, 
respectively. The MFL model yields Ff ~ constant and 
F^ ~ y/r, respectively. None of the models considered 
have presented analytic forms for the scattering rate as 
a function of doping, and presumably in all models F^ 
would be expected to be small in overdoped systems. 

It is often useful to look at the "scattering ratio" 
Ff (r)/Ff (T) ~ pc{T)/pab{T) ~ T-™. The models dis- 
cussed give TO = 1/2,3/4, and 1 for MFL, "hot", and 
"cold" spot models, respectively. These preceding expo- 
nents are summarized in Table I. As can be seen from Fig. 
([T|) , all of these models can qualitatively describe the ex- 
perimental results for overdoped systems, but important 
deviations occur for optimal and underdoped systems. 
The "hot spot" model yields a stronger temperature de- 
pendence however than that seen for the Big Raman and 
c-axis conductivity. 



appearing at a high energies which merges into the su- 
perconducting feature at Tc. It is not immediately clear 
whether this occurs in Raman data due to the nature of 
the 600 cm~i peak2^. 

In the spin and/or charge precursor scenario, 
anisotropic SDW and/or CDW fluctuations strongly af- 
fect the integrity of qp excitations near the BZ axes^^^. 
Relatedly strong electron and Umklapp scattering in par- 
ticular due to the nearness of a nesting condition can 
drive FS topological changes near the BZ axes or "hot 
spots" which preservoSi or lower— the symmetry of the 
FS. 

It is clear that the pseudogap is a manifestation of 
strong correlations regardless of which scenario is con- 
sidered. Thus we take a simple approach and relate the 
pseudogap to a correlation gap as a precursor to the 
Mott insulating phase regardless of any underlying or- 
der, characterized by the development of Gmc involving 
large energy transfers across a precursor Mott gap. The 
gapping can thus be crudely understood as the loss of 
well defined qps located near the (tt, 0) regions of the FS, 
implying that the coherent part of the Green's function 
diminishes away from the BZ diagonal. 

Therefore in what follows the role of anisotropy in the 
qp residue Z is explored in a simple effort to model the 
effect of a loss of qp transport for the "hot" qps with 
decreasing p. Taking Zi,{T) = ^^g-^^"/^) ™"'(2<[.) a 
phenomenological model of angular dependent gapping 
of qps with an energy scale Eg , the integrals are straight- 
forward and the result can be expressed analytically in 
terms of a degenerate hypergeometric function of two 
variables: 



2rc(T) j $1 
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C. Pseudogap 

The upturn of both Fi(T) and Pc{T) at low tempera- 
tures for optimal and underdoped systems is indicative of 
gapped qps and connected to an anisotropic pseudogap 
largest near the BZ axig^. A major issue^^ is whether 
the pseudogap is caused by pairing without long-range 
phase coherence or due to loss of well defined qps at the 
FS related to the formation of a precursor Mott gap, or 
spin-density and/or charge-density wave states, for ex- 
ample. 

In the former case, the superconducting gap ampli- 
tude closes at T* while strong phase fiuctuations force 
the superfluid density to appear at T^^^. In more exotic 
phases emerging from Z2 gauge theories, electrons frac- 
tionalize away from the BZ diagonals, spinous become de- 
confined and holons condense and become gapped^"^. In 
these scenarios one might expect a feature in the spectra 



For almost all temperatures, the function can be ac- 
curately described as the previous results Eqs. H12t - 
I15|l with the sole exception that a'^^{T) Eq. H13|) and 
dx'h, {T)/dn Eq. lUni are muhiplied by e^^Sg/T rpj^^^^ 
we note that if qps located near the BZ axis become 
gapped or lose their spectral weight at the Fermi level, 
the Big Raman slope and c— axis conductivity will show 
activated behavior while the i?2g Raman slope and the 
in-plane conductivity would continue to show metallic 
behavior. This is qualitatively the situation found for 
doping levels below pc ~ 0.22 in all the cuprates. 



D. Comparison with experiments 

The data for the Raman derived scattering ratio for 
Bi-2212 is shown in Fig. |21 The data are derived from 
the measurements shown in Fig. ^ The ratio derived 




T[K] 

FIG. 2: A log-log plot of the Raman derived "scattering ra- 
tios" rf/r|* (defined in the text) for Bi-2212 in Refi^ for 
underdoped (circles, m = 0.36), optimally doped (squares, 
m — 0.76), slightly overdoped (diamonds, m — 0.614) and 
appreciably overdoped (triangles, m = —0.13) samples shown 
in Fig. respectively. The exponent m is determined from a 
least-squares fit to T"™. 



from the measurements on three differently doped sam- 
ples of Y-123 are shown in Fig. E| For Bi-2212 the ratio 
slightly increases (to < 0) with temperature for appre- 
ciably overdoped systems, in agreement with the results 
obtained for La-214S. For decreasing doping p in both 
Bi-2212 and Y-123, the exponent to is positive and in- 
creases as the "hot" qps become gapped and the "cold" 
qps do not appreciably change. The large variation of the 
data from the underdoped Bi-2212 sample is due to the 
small intensity at low frequencies from which the slope is 
derived. Apart from this sample however a power-law fit 
adequately describes the data for both compounds. Near 
optimum doping both the MFL and "cold spot" model 
give reasonable agreement for the "scattering ratio" while 
on the underdoped side the "cold spot" model gives an 
exponent of 1 in agreement with the data on Y-123 and 
in rough agreement with the data on Bi-2212. An ex- 
ponential dependence on temperature has been used in 
Refsji for the resistivity ratio to determine the magnitude 
of a pseudogap, for example. We note that the Raman 
measurements are not yet of sufficient precision to deter- 
mine Eg from a fit since a straight line fit works well as 
shown in Fig. |21 The curvature may be obscured by the 
small signals measured at low frequencies however. More 
accurate data would be very useful. 

Recent ARPES measurements have revealed that the 
qp self energy may not be as anisotropic as determined 
earlier due to the more accurate detection of bilayer split- 
ting near the BZ axes^. In addition, a more quantita- 
tive investigation of the qp self energy derived from re- 
cent ARPES measurements on overdoped and optimally 




T[K] 

FIG. 3: A log-log plot of the data obtained in Ref^ for 
the Raman derived "scattering ratios" Ff/F^ (defined in 
the text) for YBa2Cu306.5 (circles, m — 1), YBa2Cu306.93 
(squares, m — 0.71), and YBa2Cu3 06.98 (diamonds, m = 
0.57), respectively. 



doped Bi-2212Sa has been used to argue that agreement 
with the magnitude and temperature dependence of in- 
plane resistivity measurements on similar compounds can 
only be obtained if the transport scattering rate has no 
contributions from Ths and is given solely by an MFT 
dependencaiii on all regions of the BZ. A similar con- 
clusion has been reached regarding the self energy and 
optical conductivity^^. It is however important to note 
that the magnitude of the derived resistivity agrees with 
experiment to within a factor less than two for temper- 
atures between 100-300K. A more or less isotropic qp 
self energy cannot be reconciled with the Raman data 
unless vertex corrections are brought into play. If the 
experiments are taken at face value, this would imply 
that vertex corrections affect the Raman response and 
the conductivity in opposite ways: Raman vertex correc- 
tions would act to enhance qp scattering near the BZ axes 
while conductivity corrections would decrease it. This 
has the important consequence from the symmetry ar- 
guments presented above that the c— axis conductivity 
should then be more coherent in contrast to experiments. 
This lack of consistency between ARPES, ab- and c- axes 
conductivity, and i?ig,2g Raman must then be traced to 
the detailed role of vertex corrections for each response. 
This work is currently in progress. 

It is important to point out that the results obtained 
on La-214 are qualitatively different from Y-123 and Bi- 
2212 in underdoped systemsiS^. For Lai.9Sro.iCu04, a 
clear Fermi liquid like peak develops at low frequencies 
in the Big channel which sharpens as temperature is low- 
ered so that dx"{T)/dil falls with decreasing tempera- 
ture, similar to the behavior of the i?2g channel in Y-123 
and Bi-2212. These features appear more or less con- 
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tinuously with doping. However, the peak in the B2g 
channel seems to mimic the response in Y-123 and 
Bi-2212. We note that this is consistent with ARPES 
in which a more smeared spectral function is seen for 
{tt/2, 7r/2) rather than (tt, 0) crossings^^. Recently strong 
far infrared peaks have been observed in ab-plane optical 
response™ in La2-3;Sr2;Cu04 for x = 0.05 — 0.19 which 
follow a dependence on x consistent with a coexistence 
of charge stripes and antiferromagnetic domainsS^. Sim- 
ilar strong far infrared peaks have also been observed in 
BisSrsCuOe (Bi-2201)'^^ and interpreted^ in terms of in- 
stabilities of a Fermi liquid to charge ordering. While this 
interpretation is still open to questions, both of these ob- 
servations can be reconciled with Raman scattering mea- 
surements if the stripes were aligned solely along the Cu- 
O bond directions. Whether the stripes are conducting 
or insulating, and whether they are static or dynami- 
cally fluctuating, the B2g Raman response would have a 
polarization component perpendicular to the stripes and 
thus would project onto incoherent qp transport channels 
while the Big would have a finite project of both the 
incident and scattering polarization light vectors along 
a sector of coherent, conducting excitations consistent 
with observations. These simple symmetry considera- 
tions would change if the stripes were thought to be fluc- 
tuating in various different orientations or rotated by 45 
degrees, as evidence suggest they might for more under- 
doped samples. More data and further calculations are 
essentially needed to clarify this point. It is an important 
and open issue to understand why this occurs for a wide 
range of doping in La-214 and not Y-123 and Bi-2212. 

We note that only limited experimental information 
exists concerning c-axis Raman measurements due to the 
surface problems, but recently Quilty et al. have shown 
that the low frequency c-axis Raman spectral weight in 
YBa2Cu408 depletes as temperatures are lowered^*^. In 
conjunction with the spectral weight depletion at low 
temperatures seen in Big measurements on the same 
compounds^, the admittedly limited experimental evi- 
dence is also consistent with Aig_zz and Big scaling. 
More data would of course be useful to check this fur- 
ther. In this regard it should be mentioned that there is 
recent evidence that the c-axis Raman may shed light on 
a Raman active c-axis plasmon^^. It would be extremely 
useful to examine whether the plasmon would violate the 
scaling relation or could possibly lead to a mode-coupling 
which pushes the plasmon into the Big channel. 



In the absence of vertex corrections, the expressions 
for the Raman response and the optical conductivity in 
the static limit are given in terms of the Nambu Green's 
functions as: 



dx''{T)/dn 



dx df{x) 

TT dx 



= 2 



k 



(20) 



G'^{Kxr±G';{Kxr} 



Go(k,w)fo 



Here G'(k, w) = — — ——^ — tttt^ 
Gi(k, a;)fi + G3(k, a;)f3 with the renormalized quanti- 
ties determined from the Pauli components of the self 
energy as Co = uj — So(k, oj), e(k) = e(k) -I- E3(k, oj), and 
A(k) = A(k) + Ei(k,(::;). 

It is well know that vertex corrections appreciably al- 
ter universal results and the Wiedemann-Franz law for 
d— wave superconductorgSii^ii^. In addition, they are 
crucially important for describing the back-flow needed 
to restore gauge-invariance in the superconducting state 
and appreciably alter the fully-symmetric Aig response 
over a wide range of frequencies'"^. Again we neglect 
them to exploit simple symmetry considerations. There- 
fore we only consider axx , f^zz and the Big and -620 Ra- 
man response. The reader is referred to RefsiSi2aiSi2& 
where these issues have been addressed at length. 

The self energy is usually broken into an inelastic 
term, such as due to phonons or spin-fluctuations, and 
an elastic term due to scattering from impurities: S = 
j^meiasuc _^ Y^eiasUcT^^ gjj^(,g ^^le integrand in Eq. ^ 

is weighted out for small frequencies and since S"3(k,a)) 
coming from inelastic scattering are odd functions of fre- 
quency while I]Q(k,ii>) is even, we only retain Sq. If 
one considers s— wave impurity scattering in the Born 
or unitary limit, then Yjf'^f*-'^'^ can be neglected as well. 
However generally in other limits and in particular if the 
impurity potential is anisotropic as it should be in corre- 
lated systems, one must keep these terms as well2*2^i2i. 

In the next subsection the role of disorder in determin- 
ing the asymptotic low temperature limit of the results 
functions is considered, and then in the following subsec- 
tion, inelastic scattering from spin fluctuations is used to 
determine the full temperature dependence below T^- 



Disorder 



III. SUPERCONDUCTING STATE 
A. Formalism 

We now consider how anisotropic transport in the nor- 
mal state may be reflected in the superconducting state. 
In particular we would like to address whether the vari- 
ant of the SS relation presented in Eq. Ijltill holds in the 
superconducting state. 



We first consider scattering from point-like impurities 
to determine the low temperature limit of the response 
functions in the superconducting state. For s— wave im- 
purity scattering a) = uj-V-^^^, with 50 = T( c5^-Aii(k) )' 
r = , rii the density of impurities, and c the phase 
shift"*^. The self energy is determined self-consistently 
for temperatures below T* rii due to the formation of 
a bound-state impurity band at the Fermi level. In this 
limit, the solution may be expanded for small frequencies 
a.s uj = auj + ijo + ihuP', with a, 6 and 70 determined from 
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the impurity concentration and magnitude of the phase 
shifts. Performing the standard integrals in Eq. I|2U|) 
yields in the limit of low temperatures T <^T* 



<,^(r«T*) 
ax" (T « r*)/ai] 



dx 



dfjx) 
dx 



' ' x' 



5/2,0 



2" 70'- ''7/2,1 



(21) 



with the functions 



72(k)A''(k) 
[72 + A2(k)]- 
«„(k)t;^(k)A^(k)\ 



[72 + A2(k)]- 



(22) 



and the constant c^''^ = 2,0 for the Raman response and 
conductivity, respectively due to the different coherence 
factors. Eqs. (^11 - 1^ reduce to those found in Refi2£ for 
the case of the a — h plane conductivity. The functions 
/jj'^ are straightforward to compute for a cylindrical FS 
and A(k) = Aocos(20). For resonant impurity scatter- 



mg (c = 0), a = 1/2,6 = 

self-consistently via 70 ~ 
then yield: 



870 



, and 7o is determined 



tFAo 



43 



21n(4Ao/7o) 



a;,(T«r*) = 



aL(r«T*) = ^2^^ 

TOTtAq Vp Aq 



mTrAo 

2 2 



1 + — ^ 



TT^ T 

T2 7o^ 



•21 



B2, 



{T<^T*) = 



2Nf 



'21 



TT ^ 
3672 



(r«r*) = ^6?^ln(4Ao/7o 



tAo An2 



■ _7r2T2 
^ 12 70^ 



(23) 
(24) 



(25) 



(26) 



where n is the 2D electron density. Eq. H23|l 
for the in-plane conductivity has been derived several 
^jjj^gg43j73j78j79^ and Eqs. H25I26|I for the Raman slope 
are identical to those found in Ref4^. The result for the 
out-of-plane conductivity for T = is also in agreement 
with the result from Ref.^^, but the temperature depen- 
dent variation has not been presented before. We note as 
in Refsi^2iiM2iSi2Si2a that both the in-plane conductivity 
and the B2g Raman slope are universal numbers for res- 
onant scattering independent of the strength of the scat- 
tering, while both the c-axis conductivity and the Big 
Raman slope depend on 70. The 70 dependence does not 
appear in the c— axis conductivity if the c— axis hopping 



is taken as a constant independent of direction around 
the FS^^'^". The temperature dependencies are positive 
for both the in-plane conductivity and the B2g slope, but 
are negative and identical for the out-of-plane conductiv- 
ity and the Big slope, giving a peak at zero T for the 
latter pair. We note that this result is in agreement with 
the rise of the c— axis conductivity recently observed in 
YBa2Cu3 6.95 at low temperaturesAZ,. 

In the limit of higher temperatures Tc :$> T :$> T* 
where the DOS does not have an impurity induced weight 
at the Fermi level and matches the DOS from the clean 
limit, the self consistency is not required for the self en- 
ergy and Eq. (|20|l can be rewritten as 

( <,/3(r*«r«Te) \ 



.2Nf I dx^&Im 



dx 

Va{'k)v/3{'k) 

7'(k) 



Re 



Q — i/T{x) 
x 



y/x^ - A(k)2 



>, (27) 



with 1/t(x) = — 2I]q(.t). This is a generalization of the 
results obtained in Refs.-i2ii44^ to the case of Raman and 
optical conductivity. We note that for d— wave supercon- 
ductors in the resonant limit, the impurity scattering rate 
depends strongly on frequency 



1/t{lu ^ 0) = 



TT^FAo 



I 



2uj ln^(4Ao/w) 



(28) 



as shown in Ref. , which yields 

\ dxl^{T* « T « T,)/dn 



4:Nf T2 
7r2 FAo 



dzz^ 



(e^ + 1)2 
X ln2(4Ao/zr)iJ'"-^'X^'^(0T), 

with the functions 

H<'^-^{x)^Rel^^^^M= 
\ v/a;2 - A(k)2 



\ Jx^ - A(k)2 / 



^^2 - A(k)2 

Performing the integrals gives for small x gives 

B^-<'-^{x) = 

W'^-^ix) = 




(29) 

(30) 
(31) 

(32) 
(33) 



The remaining integrals in Eq. H29|l can be easily per- 
formed; 

a' (r*«T«re) = 



10 
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FIG. 4: Temperature dependence of the in-plane conductiv- 
ity (Panel A) and the -B29 Raman slope (Panel B) for res- 
onant scattering and different impurity scattering strengths 
r/Ao = 0.004,0.008,0.016,0.024 and 0.04 (solid, dotted, 
short-dashed, long-dashed, and dotted-dashed lines), respec- 
tively, for Ao/Tc = 4. Here ao = nNpe^Vp/Ao and x" = 
x" l^phi are the dimensionless quantities shown. 
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The expression for the in-plane conductivity was de- 
rived in Ref.^^ but to our knowledge the other terms are 
new. We note the results for Uzz and axx in this limit 
differ from those obtain in Ref^, where a frequency 
independent scattering time was used rather than that 
of Eq. H28() . As a consequence they concluded that 
<Jxx,zz{T) (X nxx,zz{T) with nxx,zz{T) the normal-fluid 
density which decreases uniformly with temperature in 
contrast to experiments^. From that they concluded 
that the scattering time must be anisotropic. We note 
that any frequency dependence of the scattering time 
would qualitatively change this conclusion. 

The results of Eqs. ^2A) and (32-|23l imply that the 
SS relations in the normal state Eq. (|16|) hold in the 
superconducting state. The exponent of the low temper- 
ature rise as well as the sign of the correction do obey the 
general scaling relation, following simply from the inter- 
play of anisotropies of A(k) and the respective vertices. 

Again these relations Eqs. (|23l -24) and (32-|3SJ) would 
be expected to be violated for the same reasons discussed 
in the normal state in Section II. However additional con- 



FIG. 5: Temperature dependence of the out-of-plane con- 
ductivity (Panel A) and the B\g Raman slope (Panel B) for 
resonant scattering (c=0) and different impurity scattering 
strengths P/Aq = 0.004,0.008,0.016,0.024 and 0.04 (solid, 
dotted, short-dashed, long-dashed, and dotted-dashed lines), 
respectively, for Aq/Tc — 4. Here cro = TvNpe^v^/ Ao and 
X" ~ x" /^pbi are the dimensionless quantities shown. In- 
sets: low temperature rise of both ctz^ and dx" /dO, (normal- 
ized to their zero temperature values) with decreasing tem- 
perature. 



siderations should be mentioned here as well. It is well 
known that at low temperatures T <^ T* the a6-plane 
conductivity in Y-123 varies as T" with an exponent 
a < 1— and not predicted by Eq. and has 

been found generally to be non-universal in the zero T 
limit^"'^. While vertex corrections can address non uni- 
veral numbcrs^"^ and scattering away from the unitary 
limit changes a from iM*^, systematic agreement has 
not been reached at low temperatures. To address this 
discrepancies, recently Atkinson and HirschfeldSS have 
shown that a reduced a6-plane conductivity emerges at 
low temperatures when real-space variations of the order 
parameter in the neighborhood of the impurities and im- 
purity interference effects are consider in a semiclassical 
Bogolubov - de Gennes framework. These effects are not 
captured in the self consistent T— matrix approach and 
arc thus beyond the scope of the present manuscript. It 
is not immediately clear how the changes in (Jab{T) are 
manifest in other response functions considered in this 
manuscript and how the derived scaling relations are af- 
fected. Our approach should be valid at not too low 
temperatures where deviations of the conductivity from 
the unitary limit results are found. 

The response for T <C is calculated by numeri- 
cally solving Eq. [^Oland the corresponding self-consistent 
equations to determine the self energies. The results 
for axx{T),dx'^^^[T)/dn, and a,,(r), 9^^^^ are 
shown in Figures 01 and [S] respectively, for resonant scat- 
tering and different values of the impurity scattering 
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strengths F/Aq. Generally at higher temperatures T > 
T* all quantities increase rapidly with temperature, ris- 
ing as and T** for a^x, dxE^g/dO. and a^^, dxBiJdil, 
respectively. The rise of the c— axis conductivity and the 
Big Raman slope at low temperatures shown in the in- 
sets of Fig. El are generally on the order of a few percent 
for the parameters shown. This height rises for smaller 
values of F but onsets at smaller temperatures due to 
the concomitant reduction in T*. In particular, the rise 
and the onset of the c— axis conductivity low temperature 
maximum for YBa2Cu3 06.93'*'' cannot be adequately re- 
produced. There are as yet no Raman measurements to 
compare to, and thus it would be extremely useful to have 
data on a wide range of compounds and doping levels as 
well as a systematic check of impurity doping effects to 
test these results. 



C. Spin fluctuations 

The different rate of descent of the response functions 
below Tc has an interesting consequence on the con- 
ductivity peak seen in a6-plane measurements and the 
lack of peak seen in c— axis measurements. Ref4^ in- 
cluded inelastic scattering from spin fluctuations in RPA 
to reproduce the a6-peak in the conductivity observed in 
RefiSiii. In RefsiM3 it was shown that Eq. ^ for the 
conductivities for Tc ^ T ::§> T* may be reexpressed in 
terms of the normal qp density which can be generalized 
as 
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with 
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(37) 

the projected normal quasiparticle density. The average 
r is derived from the frequency-dependent t{ll}) and the 
superconducting DOS N{llj) as 



J dujN{uj){-df/duj)T{Lj) 
J dujN{uj)i~df/duj) ■ 



(38) 



Similarly one can re-express the Raman slopes in the 
same fashion: 



with 



qp 



(T) = j dul^Re 
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(40) 



the Raman projected normal qp density. For a d^^-y^ 
gap, from the results of Eqs. (28-31) the projected qp 
densities at low T vary as T for n'^pjU^p^^^ and for 
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FIG. 6: Frequency dependence of the imaginary part of the fo 
self energy Eg (k, oj, T) normalized to the hopping overlap t as 
a function of frequency and temperature for different points in 
the BZ. The solid line and dotted line are for T — O.STc for gap 
maximum k — (7r/a,0) and gap node (7r/2a, 7r/2a) point, re- 
spectively, while the dashed and dot-dashed lines correspond 
to the gap max and gap node points at Tc. The inset shows 
the zero frequency part of Eq as a function of temperature. 



n^p, n^p^^" , respectively. If the scattering time r were in- 
dependent of frequency, then Ugp gives the full tempera- 
ture dependence and thus axx , 9x'B2g I would vary lin- 
early with T and a^z , dxsig I would vary as . Ref 
used this result for a^z and argued that accurately fit 
the data for T > 40i4r, but they could not explain the rise 
at low T . However, the impurity scattering rate as well 
as the scattering due to inelastic collisions, such as spin 
fiuctuations, depend on momentum and strongly depend 
on both temperature and frequency. The latter is cru- 
cially needed in order to explain the peak in the a5-plane 
conductivity. 

j^gfg|42j43j74 Tjitilized calculations of the inelastic scat- 
tering due to spin fluctuations in the 2D Hubbard model 
in the Random Phase approximation (RPA) for U — 2t 
to describe the dc and IR conductivity and the frequency 
dependent Raman response. The lifetime calculated for 
U = 2t and Aq /Tc = 3 — 4^^ was found to give reasonable 
agreement with the transport lifetime determined from 
conductivity measurements in Y-122'*^ and gave reason- 
able agreement with the afe— plane conductivity peak^, 
a6-plane IR conductivity response^, and simultaneously 
the a6— plane IR and the Big and B2g Raman response 
in Bi-2212^^. We therefore use this approach to calculate 
the temperature dependence of the response functions for 
all temperatures below Tc- 

In RPA, the self energy Eq is given from the effective 
potential V as: 



3 U^xoiq,^ 
21- [7x0 (q, in)' 



(41) 
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where [/ is a phenomenological parameter [we choose [7 = 2t]. xo{^.^''■^) is the non-interacting spin susceptibiHty, 
I 



Xo 



( ^ \p [ "k,k+q /(£'k+q) - fjEk) "k,k+q 



1 - /(gk+q) ~ /(i^k) 1 - /(i?k+q) - /(jgk) 

iQ + i?k+q + -Ek ifi — -Ek+q — -Ek 



Here i?k — ^k + ^k ^^'^ the coherence factors are = 1 ± ^''^"g'^^'^^^''^'' ■ This yields a self energy 



dx ^ V""(q,2;) 



ttN ^ 2£;k-q 



Sk-q + X — 

with n the Bose factor. 



-Ek-qTQ + ek-qT3 + Ak-qTl 
— i?k-q + X — iuj 



n{x) + /(i?k-q)] 
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The imaginary part of the tq self energy Sf^k, lJjT) 
normalized to the hopping overlap t as a function of fre- 
quency and temperature for different points in the BZ is 
shown in Fig. Here we have used the band structure 
Ek as in (A4) in the appendix with t' /t = 0.45 and a 
filling n — 0.88, U = 2t, and a d^^-y^ energy gap Ak = 
Ao[cos(fc:ra)-cos(fcya)]/2 with Ag/t = 0.4 = ATc/t. The 
solid line and dotted line shows the frequency dependence 
of S" at a temperature T = O.STc for gap maximum 
k = (7r/a,0) and gap node (7r/2a, 7r/2a), respectively, 
while the dashed and dot-dashed lines correspond to the 
gap max and gap node points at T^- The differences 
for the gap maximum and gap node points are not too 
strong and can be adequately fit with a threshold behav- 
ior ^ {uj — 3A(k))^ plus a temperature dependent part 
which also depends on momentum. The inset shows the 
zero frequency part of Eq as a function of temperature. 
Except for low temperatures where the nodal properties 
of the interaction govern the behavior, the momentum 
dependence of the self energy is weak and can be ade- 
quately modelled by a temperature dependent ^ term 
plus a frequency dependent part ^ oj^. 

In an effort to address the temperature dependence of 
these quantities, we employ a simple parameterized fit to 
the numerical results for 1/t\^{u!,T) — —2Y,Q(k,oj,T) de- 
termined from Eq. and Fig. © and add that to the 
elastic contribution calculated in the last section. Assum- 
ing Matthiessen's law to hold in this case neglects vertex 
corrections and the joint influence of disorder on the spin 
fluctuations and vice-versa, but for weak disorder should 
be sufficient to capture the qualitative behavior of various 
quantities derived on the FS. 

The results for the four response functions derived from 
Eq. H20|) are shown in Figures [T] and |S1 Both the in-plane 
conductivity (Fig. [7|^) and the i?2g Raman slope (Fig. 
E^) possess a peak near T ~ O.BTc for F/Aq = 0.004 
which decreases in height and moves to higher tempera- 
tures for increasing impurity scattering. It is important 
to emphasize that this peak is not related to coherence 
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T/Tc 

FIG. 7: Temperature dependence of the in-plane con- 
ductivity (Panel A) and the -B29 Raman slope (Panel B) 
including inelastic spin fluctuations and resonant impu- 
rity scattering for different impurity scattering strengths 
r/Ao = 0.004,0.008,0.016,0.024 and 0.04 (solid, dotted, 
short-dashed, long-dashed, and dotted-dashed lines), respec- 
tively, for Ao/Tc = 4. Here ao = irNpe'^v^/Ao and x" = 
x" l^piii a-re the dimensionless quantities shown. 



effects and is a simple balance of fall-off of the inelastic 
scattering rate ~ and the rise of the impurity scatter- 
ing rate ~ 1/T at low temperatures. The sum of these 
scattering rates is multiplied by the normal qp density 
oc T at low temperatures, as shown in Eqs. and 
and therefore cjxx and Oxb^^/^^^ vary as 1/T^ for 
T* < T < Tc and approach a universal constant in the 
zero temperature limit. 

However, no corresponding peak is found for both the 
out-of-plane conductivity (Fig. and Big Raman 

slope (Fig. IHJ3), in agreement with experimental obser- 
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0.2 0.4 0.6 0.8 1 
T/Tc 

FIG. 8: Temperature dependence of the out-of-plane con- 
ductivity (Tz2 (Panel A) and the Big Raman slope (Panel 
B) including inelastic spin fluctuations and resonant im- 
purity scattering for different impurity scattering strengths 
r/Ao = 0.004,0.008,0.016,0.024 and 0.04 (solid, dotted, 
short-dashed, long-dashed, and dotted-dashed lines), respec- 
tively, for Ao/Tc = 4. Here ao = nNpe'^vl/Ao and x" = 
x" /Npbi are the dimensionless quantities shown. 



vations. The curves simply show a rapid fall off of both 
quantities for T < Tc and a small rise of both quantities 
which onsets at T* and reaches a zero temperature max- 
imum as shown in Fig. \S\ The main difference is due 
to the behavior of zone-axis projected qp density, which 
varies as at low temperatures, with a factor of T com- 
ing from the nodes and the additional coming from the 
matrix elements. This compensates the rise of the 

qp inelastic lifetime, and both azz and Oxb^ vary as 

for T :$> T* with fi dependent on the strength of the 
impurity scattering, and rise for T <C T*, as shown in 
Fig. (0. For example, for the parameters chosen in Fig. 

the exponent fj. for O.STc < T < 0.9Tc for the c-axis 
conductivity azz{T) varies from 2.7 to 3.4 for increasing 
impurity scattering. If the frequency dependence of the 
scattering rate were neglected, then a universal exponent 
/I = 3 would emergei^. Therefore it would be highly 
useful if further systematic checks were performed and 
Raman data were available to compare to the conductiv- 
ity and the theoretical predictions. 



IV. SUMMARY AND CONCLUSIONS 

In summary, based on symmetry arguments we have 
demonstrated how the relaxational behavior of the qp in 
the cuprates should manifest itself in the various experi- 
ments and how the results are expected to be interrelated. 
Therefore, a single framework may relate the optically- 
derived qp scattering rates to transport measurements 



to infer charge dynamics on different regions of the Bril- 
louin zone. Using forms for the interlayer hopping and qp 
self energy consistent with empirical evidence, a variant 
of the SS relation was shown to relate the zone-diagonal 
and zone- axis transport properties measured by DC con- 
ductivity and the slope of the Raman response in the nor- 
mal state, in agreement with experimental observations 
in Bi-2212 and Y-123, but not La-214. Violations of the 
derived scaling relations were discussed most pointedly 
in connection with the role of vertex corrections. 

The "scattering ratios" show power-law behavior for 
the Raman response which can be reasonably accounted 
for in several models near optimal doping. However no 
single model can adequately describe the data over the 
entire doping range, indicating that additional physics 
related to strong correlations is required^®. A presence 
of a pseudogap is discussed in simple symmetry terms, re- 
vealing that the Big Raman scattering and c-axis conduc- 
tivity are most affected in agreement with experiments. 
This is a consequence of a loss of qp coherence near the 
BZ axes. 

The data on La-214 over a wide range of doping is in- 
consistent with the simple models for qp scattering dis- 
cussed herein. A connection can be made between the in- 
plane conductivity and Raman response in light of stripe 
orientation. However more work is clearly needed to ad- 
dress this point. 

We note that a quantitative connection between the 
magnitude and temperature dependence of the qp self 
energy derived from ARPES, the in-plane and out-of- 
plane conductivity, and the Raman response can only be 
undertaken with an understanding of the role of vertex 
corrections. 

In the superconducting state, a similar SS relation is 
found which arises from the momentum dependence of 
the energy gap and conductivity and Raman matrix el- 
ements. In particular, we found that the a zero tem- 
perature peak is predicted to arise in azz and dxE^^/d^ 
without the presence of a maximum near 0.3Tc found 
for axx and Oxb^ /dO.. The results are in rough, qual- 
itative agreement with the available data for Uzz but 
the strength of the elastic scattering cannot simultane- 
ously account for in-plane and out-of-plane conductiv- 
ities. However, the simple model presented does not 
account for anisotropies in impurity scattering, known 
to arise for point-like scatterers in correlated materials, 
or impurity interference effects. In particular, it would 
useful to determine whether the approach followed by 
Atkinson and Hirschfeld*^ would remedy the agreement 
to the c— axis conductivity rise at low temperatures. Un- 
fortunately, Raman data in the superconducting state to 
further test the theory is lacking. In particular, it would 
be extremely useful to determine if the deviations from 
the derived SS relation observed in the normal state of 
La-214 carry over into the superconducting state. 

The agreement of the derived SS relations in both the 
superconducting and normal states with the available 
data on Bi-2212 and Y-123 indicate that the in-plane 
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momentum is at least partially conserved in c-axis trans- 
port over the entire doping range studied. This shows 
that in principle a comparison of Raman and transport 
could eventually contribute to the solution of the c-axis 
transport problem. 
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IjPijj y = -Px,-y = respectively. Lastly the bracket 
(• • •) notes a sum over the nearest neighbor Cu4s sites 
in the c— direction. Thus c— axis hopping is mediated 
by the Cu4s orbitals hybridizing with the bonding and 
anti-bonding pd bands consistent with LDA^^. 

After Fourier transforming, the Hamiltonian is H = 
J2k,a -^k,cr with 

-{2iipddl [aa;,k,crSx(k) - aj^,k,o-Sa(k)] -|- h.c.} 
-{2itpsSl^Ja^^k^„Sx{k) + ay^]^^^Sy{k)] + h.c.} 

-4ippS^(k)sy(k)[a^^^ay^k.<T + h.c], (A2) 



with SQ,(k) = sin(a/cQ,/2) and es(k) = — 2tggC0s{kzc). 
APPENDIX A Eq. IA2I can be diagonalized by defining "canonical 

fermions"— : 

We start by considering a four band model for the 

Cu02 plane with Cu^^ — Op^ ^ hopping amplitude ip^, _ . Sx{^)o,x,k.,a ~ Sy 
Cu4s - Op^ ^ hopping amplitude tps, Op^ - Op^ hopping fi(k) 
amplitude tpp, and c-axis Cu4s-Cu4s amplitude ts, re- .Sy{k)a^,u,a + Sx{k)ay^u,a 

spectively: /^k,cr = — « ~^/i(k) ^ ' 



- tpd J2 Psidl 

- hs X! ^5'(4,<T«n,<S,'T + h.c.) 

n,(5,fT 

(n.m) ,a 

where es,d — Eg^d ~ Ep represents the charge transfer 
energy from the oxygen p— to Cua^sM orbitals, respec- 
tively. Here Sn^o-j^n.o- creates an As,'idx^-y2 electron, 
respectively, with spin a at a copper lattice site n, while 
a„ 5,(T annihilates an electron at one of the neighbor- 
ing oxygen sites n -I- 5/2 determined by the unit vector 
5 assuming the four values, (±1,0) and (0, ±1). The 
overlap factors P have the following properties: ^'(i,o) = 



where ^(k)^ — s^(k) + s^(k). This gives anti-bonding, 
bonding bands hybridized with the Cu orbitals. This 
four-band model can be reduced to an effective one-band 
model by eliminating the (3 band and the two bands with 
high energies ~ e^^d- This is achieved by defining two 
other sets of canonical fermions and expanding in powers 
oi tpd,pd,ss/is,dF^- The single-band dispersion is approxi- 
mately given by 

e(k) = — 2i[cos(A:a;a) -I- cos(fcj,a)] -I- At' cos{kxa) cos(fcj,a) 
—2t" cos{2kxa) cos{2kya) 

—t± cos(/c^c)[cos(/cj;a) — cos(A:j,a)]^ — /i, (A4) 

with the identification to lowest order of t = tpp — 
ildl^d.t' = ~tpp/2 + tU&ts.t" = tlJl6e„ and t^ = 
tsstps/^l- This form for the interplane hopping can also 
be derived in the framework of the Hubbard model by 
projecting out the high-lying Cu 4s orbitals and the high- 
lying d~p spin triplets by solving the correlation problem 
within the unit cell and treating the intercell hopping as 
a degeneracy lifting perturbationi2*2k. 
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